Pion Imaging at the AGS 
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Source imaging analysis was performed on two-pion corre- 
lations in central Au + Au collisions at beam energies between 
2 and 8A GeV. We apply the imaging technique by Brown and 
Danielewicz, which allows a model-independent extraction of 
source functions with useful accuracy out to relative pion sep- 
arations of about 20 fm. We found that extracted source 
functions have Gaussian shapes. Values of source functions 
at zero separation are almost constant across the energy range 
under study. Imaging results are found to be consistent with 
conventional source parameters obtained from a multidimen- 
sional HBT analysis. 



I. INTRODUCTION 

Measurements of two-particle correlations offer a pow- 
erful tool for studies of spatial and temporal behaviour 
of the heavy-ion collisions. The complex nature of the 
heavy-ion reactions requires utilization of different probes 
and different analysis techniques in order to obtaine re- 
liable and complete picture of the system created in the 
collision. 

We will discuss application of imaging technique of 
Brown-Danielewicz which allows one to reconstruct the 
entire source function, phase space density and entropy 
from any like-pair correlations in a model independent 
way. We will present recent results of the two-pion imag- 
ing analysis performed by the E895 Collaboration and 
discuss natural connection between imaging and tradi- 
tional HBT. 



A. E895 Experiment and Analysis Details 

Detailed description of the experiment E895 can be 
found elsewhere Q. The Au beams of kinetic energy 
1.85, 3.9, 5.9 and 7. 9 A GeV from the AGS accelerator at 
Brookhaven were incident on a fixed Au target. The an- 
alyzed 7r~ samples come from the main E895 subsystem 
- the EOS Time Projection Chamber (TPC) |§, located 
in a dipole magnetic field of 0.75 or 1.0 Tesla. Results of 
the E895 ir~ HBT analysis were reported earlier ||, and 
in this paper, we focus on an application of imaging to 
the same datascts. 



B. Source Imaging 

In order to obtain the two-pion correlation function C2 
experimentally, the standard event-mixing technique was 
used |I| . Negative pions were detected and reconstructed 
over a substantial fraction of 4ir solid angle in the center- 
of-mass frame, and simultaneous measurement of particle 
momentum and specific ionization in the TPC gas helped 
to separate 7r~ from other negatively charge particles, 
such as electrons, K~ and antiprotons. Contamination 
from these species is estimated to be under 5%, and is 
even less at the lower beam energies and higher transverse 
momenta. Momentum resolution in the region of corre- 
lation measurements is better than 3%. Event centrality 
selection was based on the multiplicity of reconstructed 
charged particles. In the present analysis, events were 
selected with a multiplicity corresponding to the upper 
11% of the inelastic cross section for Au -I- Au collisions. 
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Only pions with pr between 100 MeV/c and 300 MeV/c 
and within ±0.35 units from midrapidity were used. An- 
other requirement was for each used 7r~ track to point 
back to the primary event vertex with a distance of clos- 
est approach (DCA) less than 2.5 cm. This cut removes 
most pions originating from weak decays of long-lived 
particles, e.g. A and K°. Monte Carlo simulations based 
on the RQMD model || indicate that decay daughters 
are present at a level that varies from 5% at 4 A GeV to 
10% at 8A GeV, and they lie preferentially at pt < 100 
MeV/c. Finally, in order to overcome effects of track 
merging, a cut on spatial separation of two tracks was 
imposed. For pairs from both "true" and "background" 
distributions, the separation between two tracks was re- 
quired to be greater than 4.5 cm over a distance of 18 cm 
in the beam direction. This cut also suppresses effects of 
track splitting Correlation functions were corrected 
for pion Coulomb interaction jj) . 
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FIG. 1. Measured two-pion correlation functions for Au + 
Au central collisions at four beam energies. 



Imaging |j|7[ was used to extract quantitative infor- 
mation from the measured correlation functions. It has 
been shown ^ that imaging allows robust reconstruction 
of the source function for systems with non-zero lifetime, 
even in the presence of strong space-momentum corre- 
lations. The main features of the method are outlined 
below; see Refs [||J7]|J for more details. The two-particle 
correlation function may be expressed as follows: 



C P (Q)-1 = J drK(Q,r)S P (r) 



(1) 



1 and $q is the relative wave- 



where K 

function of the pair. The source function, Sp(r), is the 
distribution of relative separations of emission points for 
the two particles in their center-of-mass frame. The to- 
tal momentum of the pair is denoted by P, the relative 
momentum by Q = pi — P2, and the relative separation 
by r. Since the correlation data is already Coulomb cor- 
rected, we may approximate the relative wavefunction of 
the pions with a noninteracting spin-0 boson wavefunc- 
tion: 



-3 

10 F~ 



• 2 GeV 
■ 4 GeV 

6 GeV 
8 GeV 



10 



X ' t 



EZ2 



10 




10 



8 10 12 14 16 18 20 

r, (fm) 



FIG. 2. Relative source functions extracted from the pion 
correlation data at 2, 4, 6 and SA GeV. 
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The goal of imaging is the determination of the relative 
source function (<Sp(r) in Eq. (0)), given Cp(Q). The 
problem of imaging then becomes the problem of invert- 
ing K(q, r) ju^. For pions, one might think that because 
Eq. ^] is a Fourier cosine transform it may be inverted an- 
alytically H to give the source function directly in terms 
of the correlation function: 

TABLE I. Radius parameters of Gaussian fits to the 
extracted source (Rs) functions and measured correlation 
functions(_Rc2) for different energies. 



E b (AGeV) 


2 


4 


6 


8 


Rs(fm) 


6.70±0.04 


6.35±0.03 


5.56±0.03 


5.53±0.05 


R C 2 (fm) 


6.39±0.2 


6.05±0.1 


5.51±0.15 


5.61±0.28 
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S P (r) = J d 3 Q cos (Q • r) (C P (Q) - 1) . 



(3) 



While this might work for vanishing experimental un- 
certainty, for realistic situations, this is a poor way to 
perform the inversion. When Fourier transforming, it is 
impossible to distinguish between statistical noise and 
real structure in the data fll|| . Instead, we proceed as 
in [fl3f and expand the source in a Basis Spline basis: 
s ( r ) = J2j S j B j( r )- With this > Eq. (§) becomes a ma- 
trix equation Ci = KijSj with a new kernel: 



Kij = J drK(Q t ,r)B 3 (r). 



(4) 



Imaging reduces to finding the set of source coeffi- 
cients, Sj, that minimize the \ 2 ■ Here, x 2 = Ylii^i ~ 

This set of source coef- 

ficients is Sj = Y. l [(K T {& 2 C)- 1 K)- 1 K T B] 3l (C l - 1) 
where K T is the transpose of the kernel matrix. The 
uncertainty of the source is the square-root of the di- 
agonal elements of the covariance matrix of the source, 
A 2 5 = (K T (A?C)- 1 K)-' S -. It has been shown § that 
for the case of noninteracting spin-zero bosons with Gaus- 
sian correlations, there is a natural connection between 
the imaged sources and "standard" HBT parameters. In- 
deed, the source is a Gaussian with the "standard" HBT 
parameters as radii: 
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where A is a fit parameter traditionally called the chaotic- 
ity or coherence factor in HBT analyses. Here, [R 2 ] is the 
real symmetric matrix of radius parameters 



[R 2 



r>2 r>2 p2 

p2 p2 p2 

rt os rt s n si 

p2 p2 p2 



(G) 



Eq. (||) is the most general Gaussian one may use, but 
usually one assumes a particular symmetry of the single 
particle source so that there is only one non-vanishing 
non-diagonal element R 2 ol . The asymptotic value of the 
relative source function at zero separation S(r — > 0) is re- 
lated to the inverse effective volume of particle emission, 
and has units of fm -3 : 



S(r -»■ 0) = 



A 



(7) 



As it was shown in § S(r -> 0) is an important param- 
eter needed to extract the space-averaged phase-space 
density. Figure ^ shows measured angle-averaged two- 
pion correlation functions for central Au + Au collisions 
at 2, 4, 6 and 8A GeV. Figure || shows relative source 
functions S(r) obtained by applying the imaging tech- 
nique to the measured two-pion correlation functions. 



Note that the plotted points are for representation of the 
continuous source function and hence are not statistically 
independent of each other as the source functions are ex- 
panded in Basis Splines [jl3| . Since the source covariance 
matrix is not diagonal, the coefficients of the Basis spline 
expansion are also not independent which is taken into 
account during \ 2 calculations. 
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FIG. 3. Values of inverse relative source functions at zero 
separation as a function of beam energy, obtained using imag- 
ing (solid circles) and HBT (open circles). 



TABLE II. Fit parameters for the Bertsch-Pratt HBT pa- 
rameterization of the pion correlation functions for E895 beam 
energies used for evaluation of S(r — » 0). 



E b (AGeV) 


2 


4 


6 


8 


A 


0.99T0.06 


0.74T0.03 


0.65±0.03 


0.65±0.05 


Ro (fm) 


6.22±0.26 


5.79T0.16 


5.76T0.23 


5.49±0.31 


R s (fm) 


6.28±0.20 


5.37T0.11 


5.05±0.12 


4.83±0.21 


Ri (fm) 


5.15T0.19 


5.15T0.14 


4.72±0.18 


4.64±0.24 


R 2 ol (fm) 


-2.43il.71 


0.43T1.03 


2.17±1.20 


-0.65il.85 
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With this technique, we may reconstruct the distribu- 
tion of relative pion separations with useful accuracy out 
to r ~ 20 fm. The images obtained at each of the four 
E895 beam energies are rather similar in shape, and upon 
fitting with a Gaussian function, values of % 2 per degree 
of freedom between 0.9 and 1.2 are obtained. Results 
of fits to the source function and correlation functions 
are shown in Table B. One can see that source radii ex- 
tracted via both techniques are similar, further confirm- 
ing the validity of a Gaussian source hypothesis. Figure || 
compares values related to effective volumes of pion emis- 
sion inferred from the standard Bertsch-Pratt pion HBT 
parametrization (open circles) with the effective volumes 
derived from the image source functions shown in Fig. ^ 
(solid circles). Essentially, this figure plots the inverse 
of the left- and right-hand sides of Eq. fj]. Results of the 
multidimensional Bertsch-Pratt fit to the pion correlation 
data have already been published in Ref. || and are re- 
produced in Table H It can be seen from Fig. || that the 
agreement between imaging and the HBT parametriza- 
tion is fairly good. Values of source functions at zero 
separations estimated via either technique are approxi- 
mately constant within errors across the 2 to 8A GeV 
beam energy range. 

In summary, we present measurements of one- 
dimensional correlation functions for negative pions emit- 
ted at mid-rapidity from central Au + Au collisions 
at 2, 4, 6 and 8A GeV. These correlation functions 
are analyzed using the imaging technique of Brown and 
Danielewicz. It is found that relative source functions 
S(r) have rather similar shapes and zero-separation in- 
tercepts S(r — > 0). Distributions of relative separation 
have been measured out to 20 fm, and the extracted 
source functions are approximately Gaussian. We have 
performed the first experimental check of the predicted 
connection between imaging and traditional meson inter- 
ferometry techniques and found that the two methods are 
in good agreement. This agreement paves the way for ap- 
plications of the imaging method to the interpretation of 
pair correlations among strongly interacting particle such 
as protons, antiprotons, etc. Values of source functions 
at zero separation which are related to the pion effective 
volumes of emission are almost constant across the range 
of bombarding energies under study. 
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